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Abstract. Many-particle confinement (localization) is studied for a ID system of 
spinless fermions with nearest-neighbor hopping and interaction, or equivalently, for 
an anisotropic Heisenberg spin-1/2 chain. This system is frequently used to model 
quantum computers with perpetually coupled qubits. We construct a bounded 
sequence of site energies that leads to strong single-particle confinement of all states 
on individual sites. We show that this sequence also leads to a confinement of all 
many-particle states in an infinite system for a time that scales as a high power of 
the reciprocal hopping integral. The confinement is achieved for strong interaction 
between the particles while keeping the overall bandwidth of site energies comparatively 
small. The results show viability of quantum computing with time-independent qubit 
coupling. 



PACS numbers: 03.67.Lx,72.15.Rn,75.10.Pq,73.23.-b 



Particle confinement and quantum computing 



2 



1. Introduction 

Spatial localization of excitations is central to quantum control and measurement. It 
is often required both for addressing excitations individually and for determining their 
states. Localization is particularly important for quantum computing. A quantum 
computer (QC) is usually thought of as a system of physical objects that can be modeled 
by two-level systems, qubits. The objects can be individually accessed to perform 
quantum operations. They also provide a natural physical basis for measurement. In 
many proposed implementations of QCs the qubit-qubit interaction is not turned off 
[T]- [Hj. Generally, the interaction leads to hopping of excitations between the qubits. 
The hopping is not related to relaxation, but they both may be devastating for QC 
operation. Therefore eliminating excitation hopping, or at least significantly slowing 
down the hopping rate, is a prerequisite for quantum computation. 

The problem of localization has a long history in physics, most notably in condensed 
matter physics, where localization is due to broken translational symmetry. Here it 
is usually discussed for particles. The relation to the case of qubit excitations is 
particularly simple for one-dimensional systems. The presence of a particle on site 
n corresponds to the nth qubit being excited, and the particle energy is the excitation 
energy. Much work has been done on single-particle localization. One of the best known 
examples is Anderson localization due to random disorder jOj , which has been studied in 
depth and by now is well understood. It has been demonstrated also, both theoretically 
and experimentally JB], that single-particle localization occurs also in a different 
type of systems, where the symmetry breaking is due to an incommensurate periodic 
potential. 

The work on localization in systems with many excitations is more recent [H^ . It has 
demonstrated that the particle-particle interaction can significantly affect localization 
We note that numerical studies of this problem are not very helpful, because 
the Hilbert space of a many-particle system is exponentially large, with the exponent 
proportional to the number of particles, and therefore only systems with a small number 
of particles can be simulated with classical computers. 

In the present paper we are interested in a specific formulation of the problem 
of localization. We study a many-particle (or equivalently, many-excitation) system 
and consider localization of all states in such system. More than that, bearing in 
mind applications to quantum control and quantum computing, we study strong on- 
site localization, where the wave function of each particle is mostly confined to one 
site. This is a stronger requirement than exponential decay of the wave functions at 
large distances. Formally it corresponds to the localization length being less than the 
intersite distance, for all states. We call such localization on-site confinement or simply 
confinement, for brevity. 

Many-particle confinement does not arise in a disordered system with bounded 
random site energies jTHI. Indeed, consider a chain where particles occupy sites. For 
short-range hopping an A^-particle state is directly coupled to ~ A^ other A^-particle 
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states. With probability oc one of them will be in resonance with the initial state, 
provided the site energies are uniformly distributed over a finite-width band. This 
means that the energy difference between the two states will be less or of the order of 
the intersite hopping integral J (we set h = 1). Resonating states are hybridized over 
time ~ J~^. 

In contrast to random disorder, in a QC the disorder can often be constructed 
at will, because site excitation energies can be individually controlled. This makes 
it possible to look at many-particle localization from a new perspective. Not only is 
"localization by construction" necessary for quantum computing, but in itself it is an 
exciting area of applications of QCs, which is of special interest to condensed-matter 
physics. This is because the very possibility of localization of all stationary many-particle 
states is an open question, as the distance between the energy levels is exponentially 
small for a large number of particles. 

In what follows we extend the previous work [17^ to show that, for an infinite 
ID system, many-particle confinement can be obtained by constructing a sequence of 
site energies. The proposed sequence has a comparatively narrow bandwidth. This 
is physically reasonable and is also important for applications in quantum computing, 
because the qubit tuning range is limited. Additionally, in a QC a smaller bandwidth 
is desirable for a higher speed of quantum gate operations. 

Rather than studying the stationary states, we put emphasis on the localization 
lifetime, i.e., the time during which all many-particle states remain confined to the 
initially occupied sites. We show that this lifetime can largely exceed J~^. We also find 
that, for moderately long sections of the constructed energy sequence, all stationary 
many-particle states can be strongly localized. 

The idea of our approach is to map the problem of the lifetime of localized states 
on the problem of many-particle scattering. We first construct site energies that lead 
to strong confinement of all single-particle states. Then delocalization of many-particle 
states corresponds to resonant scattering in which, as a result of the interaction, particles 
move away from the initially occupied sites. For real scattering to happen, the final 
configuration must be close in energy to the initial configuration. The lifetime of a state 
is proportional to the inverse of the matrix element of hopping onto a resonating state. 

A large localization lifetime is obtained if resonant scattering between many-particle 
configurations requires several virtual single-particle transitions via nonresonant states. 
The amplitudes of such transitions are small. The lifetime is determined by their 
minimal number, which we call the "order" of the transition. When it is limited to 4, it 
is sufficient to consider transitions that involve no more than 4 particles. The approach 
developed in this paper allows us to obtain long lifetime not for comparatively weak 
PTj but rather for strong interparticle coupling. As we show, this can be done without 
a significant increase of the bandwidth of site energies compared to the weak-coupling 
case. The strong-coupling confinement is of interest, in particular, for some models 
of QCs, like a QC based on electrons on hehum, where the dimensionless coupling 
parameter A may be as large as ~ 10 [18\ . 
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The paper is organized as follows. In Section II a qualitative picture of localization 
is provided. Section III contains a description of the effective localizing sequence J7] 
and a discussion of single-particle localization. In Sec. IV the localization lifetime 
is introduced and many-particle scattering is described. In Section V the bands of 
combination energies for many-particle transitions are analyzed and the conditions for 
eliminating resonances for transitions up to 4th order are found. Section VI contains a 
brief discussion of the results. 

2. Localization by constructed disorder: preliminaries 

Finding the ways of localizing excitations has been appreciated as an important problem 
of quantum computing. A natural formulation here is to map a QC onto a system 
of interacting spins, with S = 1/2. The spins are in a magnetic field, and the 
qubit excitation energy is the spin Zeeman energy. Several approaches to solving 
the localization problem have been put forward. An interesting idea substantiated by 
extensive numerical calculations is to arrange the Zeeman (site) energies into a ladder, 
so that none of them is in resonance with each other |T^. The bandwidth of site 
energies in this case linearly increases with the size of the system, and localization was 
studied for the Ising interaction between the spins. Another interesting idea is to create 
interaction- free subspaces pUj. The qubits that form these subspaces are separated 
from each other by the isolators, which are qubits prepared in a special state. The 
bandwidth of site energies remains finite in this case even for an infinite system, and the 
spin-spin interaction is not limited to the Ising one. However, the preparation requires 
that initially the interaction between the qubits is turned off. 

We propose a different approach, largely inspired by the results on Anderson 
localization in condensed- matter physics. We show that, by constructing site disorder, 
confinement can be achieved for a finite energy bandwidth and for the nearest-neighbor 
interaction that has a general form and is never turned off. 

Constructing a sequence of site energies e„ that leads to single-particle (single- 
excitation) confinement is quite straightforward. If all e„ are the same, as in Fig. ^^a), 
intersite hopping leads to formation of a band of fully delocalized states, with energy 
bandwidth proportional to the hopping integral J. 

The first step towards localization is to detune neighboring sites from each other 
by an energy h ^ J, see Fig. ^b). In this case nearest neighbor hopping becomes 
nonresonant. However, this is insufficient for localization, even though in quantum 
computing hterature there are claims to the opposite. Indeed, a particle can make a 
virtual transition to a neighboring site and from there to the resonating next nearest 
site. The corresponding transition matrix element is J'^/h. Therefore the states will be 
again band-like, with two bands of width oc J'^/h, which are centered at ±h/2. 

A natural next step is to detune the energies of 2nd neighbors. The detuning in 
this case should again exceed the matrix element of effective hopping. However, since 
the latter is ~ J^/h ^ J, the detuning can be smaller than h hj a factor ~ J/h, see 
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Figure 1. The successive detuning of site energies £„ in order to achieve single-particle 
localization in a chain, (a) All energies are the same, e„ = 0. (b) Nearest neighbors 
are detuned from each other by an energy h. (c) Second, fourth, and higher-order 
neighbors are detuned from each other, with the detuning that becomes smaller and 
smaller with the increasing inter-site distance. The data in (b) and (c) correspond to 
the energy sequence (0) with a — and 0.1, respectively. 

Fig. life) . At the next step it is necessary to detune 4th neighbors, but here the detuning 
can be smaller by an extra factor J/h. The process should be continued. It leads to 
a geometric-series type level detuning, with a finite overall bandwidth of the energy 
spectrum. A specific example is discussed below. 

A significantly more complicated situation arises for a many-particle system. Here, 
one should tune away from each other not only site energies, but also their combinations. 
This is necessary in order to eliminate resonant multi-particle transitions. The total 
number of states for particles on an L-site long chain is L!/A^!(L — A^)! oc 4^ for A^ = 
L/2, and therefore for large A^ the distance between the energy levels is exponentially 
small. This makes the problem of eliminating resonances hard. Nevertheless, as 
mentioned above, one can eliminate resonances that require at least a certain minimal 
number of virtual transitions to nonresonant states, leading to a long time over which 
all many-particle configurations remain confined to their sites. 

3. Localization of stationary states 

3.1. One-parameter sequence of site energies 

A ID spin system or a qubit chain can be mapped onto a system of spinless fermions 
via the Jordan- Wigner transformation [21]. The Hamiltonian of the fermions is 

H = Ho + H,, (1) 
Hi = JA \ ^ al^al^.ian+idn- 

< 'n 

Here, the fermion creation and annihilation operators. The site fermion 

energies e„ are the Zeeman energies of the spins or the excitation energies of the qubits. 

The spin-spin interaction leads to two effects, in terms of fermions. One is the 
intersite fermion hopping, with the hopping integral J. It is determined by the XY spin- 
spin interaction S^S^^-^^ + SyS^_^_l or, in the case of qubits, by the interaction between the 
qubit transition dipoles. Fermion hopping is a single-particle effect: the XY spin-spin 
interaction does not lead to interaction between the fermions. A different role is played 
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by the Ising spin-spin interaction S^Sl^^^, or, in the case of qubits, by the change of 
the transition energy of a qubit depending on whether or not the neighboring qubit is 
excited. It leads to the interaction of fermions, which is described by the Hamiltonian 
Hi. The parameter A characterizes the anisotropy of the spin-spin interaction: for the 
isotropic Heisenberg model A = 1. For concreteness we set J, A > 0. 

A simple bounded sequence of site energies that implements the construction 
described in Sec. 2 for a semi-infinite chain with n > 1 has the form |17j 



h 

2 



k=2 



(2) 



where [-J is the integer part. 

Besides the energy scale h, sequence (j21) is characterized by one dimensionless 
parameter a < 1. The overall energy bandwidth is W < h/{l — a). For a -C 1 
the energies El) form two subbands centered at ±h/2. Numerically, the subbands remain 
well separated for a < 0.4. 



3.2. Single-particle localization 

The expression Q has a special symmetry: as n increases, the coefficients at any given 
power a'^ are repeated with period 2(g + 1). This symmetry was used to obtain an 
analytical expression for the amplitude i^„(m) of a single-particle transition from site 
n to site n + m iI7j. The amplitude Kn{rn) quasi-exponentially decays for all states in 
a semi-infinite chain. For small a 



Kn{m) = a-""' (J/2/i)™. 
The exponent u is determined by lim^ 



(3) 



^ log[Kn{m)]. It is not the reciprocal 
average decay length: it varies from site to site, u = z^(n), as shown in Fig. |21 




300 



Figure 2. The exponent ly = v{n) of the a dependence of the transition amphtude 
Kn{m) lO as a function of site n. The data refer to the distances m = 200 (diamonds) 
and m = —200 (circles). The dashed hnes show the analytical limits on v for m —>■ oo. 
The values of v for the same site but for positive and negative m are different but lie 
within the same bounds. 



The values of u are bounded to a narrow region centered at = 
0.89 < u < 1.19. For estimates one can use = 1, i.e., 

Knim) ^ fsT™, K = J/2ah. 



1, with 
(4) 
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Along with decay "forward" , i.e., for m > [T7j, we have studied decay "backward" , 
i.e., for m < 0. Because of the symmetry of the sequence we expect that the decay 
backward should be also quasi-exponential in m, and the bounds on u should be the 
same as for positive m. The numerical data in Fig. |21 substantiate this conjecture. 
Therefore in the estimates given below we use Eq. Q in the form Kn{m) ^ i^l™-!, which 
applies to both positive and negative m. 

The characteristic decay length of the transition amplitude is 1/| Ini^'l. For K <^ 1, 
i.e. J ^ 2ah, all single-particle states are strongly localized and essentially confined to 
single sites, with small tails on neighboring sites. 

3.3. Inverse participation ratio 

The extent to which stationary states in a system are confined can be conveniently 

characterized by the inverse participation ratio (IPR). For an A^-particle eigenstate 
li'Nx) (A enumerates the eigenstates) it is given by 

^^^= (E . . \{^n,...n^\^/^Nx)\'y\ (5) 
V'^ — ^ni<...<n]v / 

where |$ni...ra]v) is an on-site A^-particle wave function (quantum register) with particles 
on sites rii, . . . , n^. 

If the stationary state lipNx) is fully localized, all terms in the sum (0) are equal 
to except one, which is equal to 1, and Inx = 1- In the opposite case of delocalized 
eigenstates all matrix elements in the sum (jSj) are small, and I^x ^ 1 for a large system. 
The on-site confinement that we are interested in corresponds to Jata being close to 1 
for all states A. 




Figure 3. The maximum inverse participation ratio /jv max for one particle {N = 1) 
in an L-site chain. The energies e„ are given by Eq. and h/ J — 20. The data refer 
to the segments of the chain with 1 < n < L. The plateau and sharp isolated peaks in 
the region 0.03 < a < 0.1 result primarily from pairwise hybridization of states. An 
example of the peak for L — 300 is shown in the inset on the linear in a scale. In the 
whole region 0.1 < a < 0.4 the IPR Ii max ~ 1.003, which indicates strong confinement 
of all stationary states. The asymptotic localization threshold a — J/2h is shown by 
the dashed line. 
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In Fig. 121 we show the results of the numerical calculation of the maximal IPR 
In max = maxA/ATA for energy sequence (j21) for one particle, = 1. In agreement with 
Eqs. dni), 0, all single-particle states are strongly localized for a > 0.1, where K < 0.25. 
It is seen from Fig. El that Ijy max sharply drops down for a ~ 0.03, which is close to the 
asymptotic value J/2h = 0.025 given by the condition K = 1. For 0.03 < a < 0.05 the 
IPR is close to 2 for all L, which indicates that only a few on-site states are appreciably 
hybridized; for example, Iix = 2 if two on-site states are hybridized with equal weight 
in the wave function 

Depending on the chain length, pairwise hybridization also occasionally occurs for 
particular values of a in the range 0.05 ^ a < 0.1. An example is the peak for L = 300 
shown in the inset of Fig. El It is due to hybridization of the states on sites 296 and 
300 Tf|. The hybridization happens because the interaction Ja|^a„+i shifts the site 
energies by ~ J'^/h. The shift is different for the state on the boundary, since it has 
less neighbors. In the particular case of the chain with L = 300 a simple calculation 
shows that this difference leads to resonance between the aforementioned on-site states. 
For the chosen J/2h in the region a < 0.06 the wave functions have longer tails than 
what follows from the asymptotic expression (jHl), and therefore the states away from the 
boundary are still affected by it. This leads to the peak of In max for L = 500, a = 0.056 
seen in Fig. El Changing L in the region a < 0.06 leads to disappearance of some peaks 
and onset of other peaks, with hybridization of different states. It is important, however, 
that for 0.4 ^ a ^ 0.1 these peaks disappear and all states are strongly localized. 

We have also studied ^\ a chain of length L = 12 with N = L/2 excitations (where 
the number of states is maximal, for a given L). For the section of the chain 1 < n < L, 
the maximal IPR /g max < 1-02 for h/J = 20, in a broad range 0.2 ^ a < 0.4, except 
for sharp isolated peaks for certain a. This indicates strong confinement of all many- 
particle states. However, the IPR for large L and could not be studied because of 
the large number of states. 

4. The localization lifetime 

Confining a many-particle system is much more complicated than a single-particle one. 
It requires eliminating all possible many-particle combinational resonances 6^ -|- . . . + 

~ ^mi + . . . + Bruk- Therefore we approach the problem of strong many-particle 
confinement from a different perspective. Rather than studying stationary states, we 
will consider the time over which interacting particles leave the sites where they were 
placed initially. We call this time the localization lifetime, tioc- 

The decay of on-site states is not exponential, at least for not too long times, and the 
time tioc is not the decay time. On-site states evolve primarily due to resonant transitions 
to other on-site states with the same energy. The transition rates between different 
resonating many-particle states vary strongly provided the single-particle states are all 
strongly localized. First there occurs hybridization of a given on-site many-particle state 
with the resonating state that is most strongly coupled to it. Then these states may 
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further hybridize with other resonating states, but this may take much longer. For our 
purpose, the locahzation hfetime is the time of the first hybridization. In other words, 
tioc is the reciprocal maximal rate of a transition to a resonating on-site state. 

An important characteristic of the dynamics of quantum systems is the coherence 
time ~ tcoh- In any real system tcoh is finite due to thermal fluctuations, decay into 
excitations of the medium, and fluctuations induced by external noise. We are interested 
in strong coherent confinement, in which case tioc calculated in the neglect of decoherence 
is ^ tcoh- We note that, in fact, decoherence can increase the time particles spend on 
their sites by suppressing coherent quantum transitions; however, this suppression is 
outside the scope of the present paper. 

In systems proposed for quantum computing J and J A characterize the rate of two- 
qubit operations, and typically the number of operations that can be performed over 
the coherence time is < 10^. Therefore, for practical purposes, in a quantum computer 
excitations will be locahzed provided tioc exceeds J^^ or (JA)^^ by a factor 10^. This 
condition must be satisfied for all many-excitation states and should apply to an infinite 
system for the quantum computer to be scalable. 

4.I. Many-particle scattering 

Even though the intersite transitions themselves are single-particle, resonant many- 
particle transitions may happen even where single-particle resonances are eliminated. 
This is a consequence of the energy exchange between interacting particles and also 
the energy change due to the interaction. The interaction (0) is nearest-neighbor two- 
particle. Therefore a many-particle transition results from a sequence of single-particle 
intersite transitions and scattering of particles on neighboring sites. In what follows we 
assume that, even though the particle-particle coupling may be strong, A 3> 1, still 
A <^ h/J. Then nonresonant many-particle effects can be considered by perturbation 
theory. 

A simple example of a Feynman diagram that describes a two-particle transition 
is shown in Fig. |3] In this example particles located initially on sites (n, n + 3) first 
make two single-particle transitions {n, n + 3) {n,n + 2) —>■ {n,n + l), then experience 
scattering by each other and exchange their energies, and then make two more single- 
particle transitions {n,n + 1) —>■ {n,n + 2) —>■ [n — l,n + 2). The overall sequence of 
transitions will be resonant if e„ + Sn+s ~ £n-i + £n+2- 

An alternative approach to the analysis of many-particle scattering is based on 
diagonalizing the single-particle part of the Hamiltonian ^\ . Since all single-particle 
states are strongly confined, the resulting states are very close to on-site states. However, 
the particle-particle interaction acquires a small off-diagonal long-range part. It is clear 
in this approach that it is the interaction that leads to resonant many-particle scattering. 

^From inspection of the diagram in Fig. |3]it is obvious that, for a many-particle 
transition to occur, the particles have to be brought to neighboring sites. This requires 
single-particle transitions via nonresonant sites. For a two-particle resonant transition 
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Figure 4. A diagram for a two-particle transition from sites {n,n + 3) to (n — 1, n + 2). 
The crosses indicate single-particle transitions to neighboring sites, and the dashed line 
indicates scattering of two particles located on neighboring sites. 

(^4, ^3) ^ (^1, ^2) (we choose for concreteness ks > ^4 and k2 > ki) the minimal number 
of the needed single-particle transitions is given by the parameter 

X = min(|A;i — p\ + \k2 — p — l \ + {k^ — p — 1\ + \k4^ — p\) . (6) 
p 

This is easy to see: the particles move from the initially occupied sites to neighboring 
sites {p,p + 1) and from there to the final sites; the minimum over p in Eq. (jH} is taken 
so as to minimize the total distance involved in the transition. 

For large x, the overall amplitude of the involved single-particle transitions is K'^, 
according to Eq. (j3)). Since the two-particle interaction energy is J A, the localization 
lifetime of the state {k^, k^) with respect to a transition to the state (fci, /C2) as given by 
the reciprocal transition matrix element is 

tioc ^ (JAA'-)-\ x>l. (7) 

This estimate is made assuming that the transition involves only one particle-particle 
collision. The role of multiple collisions will be studied in a separate paper [22] • 

We are interested in the localization lifetime with respect to all resonant transitions 
and for all states. Then, from Eq. (0), x should exceed a certain minimal value >imin 
for all resonant transitions. This means that all resonances separated by less than Xmin 
steps must be eliminated. j 

In what follows we will show explicitly how, in an infinite system, all resonances 
with X < 4 can be eliminated for moderately large h/J <^ 10^ even for A as large as 
A ~ 10. Then, with K ~ 0.1, the single-collision estimate ((Zj) leads to tioc > lO^(JA)^^. 

The important simplification from considering resonances with x < 4 is that the 
analysis can be limited to transitions in which only two, three, and four particles change 
their sites. We leave the analysis of three- and four-particle transitions for a separate 
paper [221; the present paper we will only give its result and concentrate on eliminating 
resonances related to two-particle transitions. 

The change of the on-site energy in a two-particle transition (k^jk^) ^ (/^i, ^^2) is 

= \£ki + £k2 ~ ^fca — £k4\- (8) 

To avoid resonances, this energy difference should largely exceed the transition matrix 
element. However, this is not enough, because of the change of the interaction energy 
that may accompany the transition. To the lowest order in J/ h it is determined by the 
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change of the number of nearest neighbors of the hopping particles. The analysis of 
diagrams of the type shown in Fig. |3 shows that, for x < 4, this number may remain 
unchanged or change by one or two, in which case the interaction energy will change by 
J A or 2JA. 

In the previous work [d we studied the case A < 1. In this case it was sufficient to 
open a gap at Se = with width > 2 J A. However, numerically the gap was relatively 
narrow, it was limited to 65/ h ~ 0.01. Therefore a direct application of the results to 
the case of large A would require using very large h/J, which is undesirable. In the 
next Section we develop an alternative approach that applies in the strong-interaction 
case. 

5. The spectrum of two-particle transition energies 

The overall structure of the differences of site energies for two-particle transitions 6e is 
fairly simple. If a = in Eq. Q, we have e„ = ±h/2. Then from Eq. (jHJ 6e takes on 
three values, 66 = 0, h, 2h. With increasing a each of these values splits into a band. 
For small a the bandwidth is oc ah and the bands are well separated. 

On the other hand, the change of the interaction energy in a two-particle transition 
is close either to zero or J A or 2 J A, as explained above. To avoid resonances related 
to this energy change, it should lie within the gaps of 6e. 

A natural way to obtain strong confinement while keeping the bandwidth h of site 
energies moderately large is the following. All transitions where the interaction energy 
change is close to zero should have 6e that is larger than this change but may be still 
much smaller than h, see below. On the other hand, all transitions where the interaction 
energy change is JA or 2JA should have 6e in the gap centered at h/2. Then the ratio 
JA/h does not have to be extremely small. Even though this is not the only option, 
we found that it allows using not too small values of a, which is advantageous in terms 
of robustness with respect to errors in the site energies. The inequality JA/h < 1/2 
imposes a much less restrictive condition on the bandwidth h than the requirement that 
2JA is smaller than the gap of Se at zero energy jl7j . 

We denote the widths of the gaps of 6e at energies 0, JA, and 2JA by 6q, 
61, and 62, respectively. These widths should exceed the matrix elements of the 
appropriate transitions. They should also be larger than the energy renormalization due 
to (possible) occupation of the second neighbors of the sites involved in a transition. 
The renormalization can be obtained by noticing that the energy change is due to 
making a virtual transition on a neighboring site, interacting with the particle on 
the next neighboring site, and making a transition back. This gives the energy shift 
~ JA{J/2hy. Therefore we need 

5o,i,2 » J'A/{2hY. (9) 
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5.1. Zero-energy gap 



Meeting the condition (jHl) for the zero-energy gap 5o requires a modification of the site 
energy sequence The most "dangerous", for x < 4, is the resonance in 5e associated 
with a transition (n, n + 1) ^ (n — l,n + 2) 17 . For this transition x = 2, whereas 
the energy difference is 5e ~ a^h with > 4 for all n = 6A; — 1. This is a very small 
5e for X = 2, because the transition matrix element is comparatively large in this case. 
In addition we have 5e ~ a^'^h when n and n + 2 are prime numbers, which leads to 
^0 — ^ in an infinite chain. 

A simple way to eliminate this resonance and to open a gap 5q is to modify energy 
sequence ^ by shifting the energy on each 6th site, 

e"^"^ = en+ {h/2)a' for n = Qk, (10) 

while keeping it unchanged on other sites, e™*^ = for n ^ Qk The parameter a' 
should be chosen in such a way as to avoid creating new resonances. 

All X < 4 transitions with small 5e can be found by modifying the diagram in 
Fig. m The rule is simply that the numbers of even and odd sites should be the same 
at the input and output [221 • The result coincides with what was obtained earlier [T7] 
using a different technique. The zero-energy gap for the sequence (fTnjl is 5o> a^h for 
a ^ a' ^ a^. This gap is illustrated in Fig. El where we plot the change in modified 
site energies 
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Figure 5. The low-energy parts of the two-particle energy differences 5e™'^ 
for the transitions with >r < 4 in which one of the involved particles is on the nth 
site (n > 2). The left and right panels refer to the energy sequence H1U|I with 
a = 0.25,0;' = 0.22 and a — 0.1, a' 1/15, respectively. In agreement with the 
analytical estimates, the zero-energy gap So/h a^. 

With the appropriately chosen a', the condition on a imposed by Eq. © takes the 

form 

> 2{J/2h)^A. (12) 

It can be seen that, when Eq. holds, not only is the many-particle renormalization 
of site energies small, but also the matrix elements of x < 4 transitions turn out to be 
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smaller than 6q. For given a and J/h, the condition (jl2j) sets the upper limit on the 
dimensionless parameter of the particle-particle interaction A. 

We have checked numerically that the gap 6q persists for a much longer chain than 
shown in Fig. 13 with n from 2 to 1000. This is sufficient to prove that the results 
apply to an infinite chain, cf. Ref. ^Zj. Indeed, in e™'^ the terms oc a' are repeated 
periodically with period 6, whereas the terms oc a"^ with different q are repeated with 
period 2{q + 1). Therefore all combinations of terms q;°, a^, . . . , a'' with g < 6 as well 
as the terms a' are repeated periodically with period equal to twice the least common 
multiple of all g + 1 < 7, which is 2 x (3 x 4 x 5 x 7) = 840. For a = 0.25 and a = 0.1 we 
have ~ 2 X 10~^ and = 10~^, respectively. The contribution to s^'^ /h of all terms 
of higher degree in a is much smaller than the gap 5q. This means that the results on 
the gap for an infinite chain will coincide with the results for a chain of 840 sites. 

5.2. High- energy gap 

We now discuss the gap in the two-particle transition energies ^e™*^ centered at h/2. 
We will be interested in the gap for all n and for all transitions with x < 4. The 
gap is determined by the widths of the energy bands of 5e^^ centered at ^e™'^ = 
and 5e^^ = h for the energy sequence (jSj), (|TO|l . Because the sequence is regular, the 
single-energy distribution has no tails and the bands have sharp boundaries. 

It is seen from Eqs. (j2I), (fTUj) that, in a two-particle transition (/C4, /^s) {ki,k2), 
the terms in , . . . , that are linear in a may all add up or pairwise compensate each 
other. The sum of all of them can be equal therefore to ±2ah, ±ah, or 0. Similarly, 
the term a'h/2 can be added or subtracted. Note that, for x < 4, only one term oc a' 
can contribute to 5e^'^. As a result, the width of the band of 5e^'^ at Se'^'^ = is 
h{2a + Oi' /2), whereas the width of the band at fe™'^ = /i is h{Aa + a'). 

The changes of the site energies in 3-particle transitions with x < 4 

I md 1 md i md md md c-™*^l 
l^fci + fcfca -r fcfca tfeg | 

also form bands centered at 0, h and 2h. To leading order in a, a', the bands at energies 
and h are bounded by {2a + a')h and h{l ± 3a ± a'). Interestingly, four-particle 
transitions with x < 4 do not change these bounds [22] • 

/^From this analysis, the gap in the changes of site energies of all x < 4 transitions 
(5£^)j<<4 is given by the condition 

2a + a' < /i-^(fe)^<4 < 1 -3a-a'. (13) 

Higher-order terms in a can be taken into account by replacing a a /{I — a) in 
Eq. fjl3|) . A direct analysis of the terms ~ shows that this gives a stronger inequality 
than what actually follows from Eqs. (j2I), (fTUI) . 

The conditions for eliminating resonant transitions with x < 4 in an infinite chain 
are given by Eq. (fT^ in which h~^{6e)^<4^ is replaced by J A and 2 J A. The distance from 
J A, 2JA to the band edges should significantly exceed the transition matrix elements 
and the energy renormalization JA{J/2hy; these quantities are small compared to ah. 
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Figure 6. The maximal inverse participation ratio /g max for a system of 6 particles 
on a 12-site chain. The site energies e™'^ are given by Eqs. (|10() with a' = 1/15. 
The left and right panels corresponds to the sections of the chain with 415 < n < 426 
and 600 < n < 611, respectively. The coupling parameter is A = 9, and h/J — 30. A 
broad region of a centered at a = 0.1 where the IPR is extremely close to 1 is seen in 
all sections of the chain that were tested ■ 

The inequalities (fT^ . (fT^ show that the range of the parameters h/J, A, and a 
where two-particle resonances with x < 4 are eliminated is fairly large. One may expect 
that, in this range, stationary many-particle states will be mostly confined to their sites, 
at least for a not too long section of the chain. This can be probed by calculating the 
many-particle IPR for a strongly coupled system. The results of these calculations for 
modified energy sequence (jlUp are shown in Fig. It is seen from this figure that, in a 
broad range of a where the localization lifetime is expected from Eqs. ( (|T2|) . (|T3l to be 
long, the IPR is extremely close to one. This shows that all stationary many-particle 
states are strongly confined. 

6. Conclusions 

In this paper we studied strong single- and many-particle confinement. We have shown 
that the previously proposed [T2j energy sequence leads to strong confinement of 
all single-particle states. The decay of the wave functions is quasi-exponential in both 
directions in the chain, with the same bounds on the decay length. 

We have also considered the problem of many-particle confinement. Here we studied 
the localization lifetime tioc during which all particles in an infinite chain remain on their 
sites. Of central interest was an extension of the results on long localization lifetime J7] 
to the case of strong particle-particle coupling, A ^ 1. We map the problem of many- 
particle transitions onto the problem of resonant scattering. A natural characteristic 
in the analysis of scattering is the difference between the initial and final site energies 
6e. Large tioc compared to the reciprocal hopping integral J^^ was obtained earlier [T7] 
when the zero-energy gap 6q in 6e for all transitions with up to 4 single-particle moves 
exceeded twice the coupling energy 2 J A. For A ^ 1 this would require a very large 
width of the band of site energies. 

We have shown that, for large A, a better option may be to decrease the bandwidths 
of the site energies of second neighbors in the chain, which leads to separation and 
narrowing of bands of 6e. Then, if J A and 2 JA are within the interband gaps of 6e, 
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resonant transitions are eliminated. As a result, the lifetime tioc scales as a high power of 
the reciprocal hopping integral, t\oc oc {JA)^^{2ah/J)^ for transitions that involve one 
particle-particle collision. The overall bandwidth of site energies required for obtaining 
a long localization lifetime is parametrically smaller than the one given by the condition 
60 > 2JA. The major limitation on the bandwidth is that it should substantially exceed 
the coupling energy. The results bear on many-electron conductivity in condensed- 
matter systems and on quantum computing with perpetually coupled qubits. They 
demonstrate the possibility to avoid delocalization of excitations in a QC without 
refocusing. 
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